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Abstract 

We report on a calculation of the full one-loop weak corrections through the order OgOw 
to parton-parton scattering in all possible channels at the Relativistic Heavy Ion Collider 
(RHIC) running with polarised pp beams (RHIC-Spin). This study extends the analysis 
previously carried out for the case of 2 — > 2 subprocesses with two external gluons, by 
including all possible four-quark modes with and without an external gluon. The addi- 
tional contributions due to the new four-quarks processes are extremely large, of order 
50 to 100% (of either sign), not only in the case of parity- violating beam asymmetries 
but also for the parity-conserving ones and (although to a more limited extent) the total 
cross section. Such ©(ogOw) effects on the CP-violating observables would be an as- 
tounding 5 times larger for the case of the LHC with polarised beams - which has been 
discussed as one of the possible upgrades of the CERN machine - whereas they would 
be much reduced for the case of the CP-conserving ones as well as the cross section. 
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The purely weak component of Electro- Weak (EW) interactions is responsible for inducing 
parity- violating effects in jet observables, detectable through asymmetries in the cross section, 
which are often regarded as an indication of physics beyond the Standard Model (SM) [1]. 
These effects are further enhanced if polarisation of the incoming beams is exploited, like at 
the BNL machine mentioned in the abstract [2, 3]. There have also been some discussions 
[4, 5] on the idea of polarising the proton beams at the Large Hadron Collider (LHC) as one 
of the possible upgrades of the CERN machine, though no proposal has been put forward 
yet. At either machine, comparison of theoretical predictions involving parity- violation with 
experimental data can be used as a powerful tool for confirming or disproving the existence 
of some beyond the SM scenarios, such as those involving right-handed weak currents [6], 
contact interactions [7] and/or new massive gauge bosons [8, 9, 10]. 

In view of all this, it becomes of crucial importance to assess the quantitative relevance of 
weak effects entering via 0{ala^) the fifteen possible 2 — > 2 partonic subprocesses responsible 



for jet production in hadronic collisions^, namely: 

99 ^ QQ (1) 

qq gg (2) 

Q9 Q9 (3) 

qg qg (4) 

qq qq (5) 

qq qq (6) 

qQ ~^ qQ (same generation) (7) 

qQ ~^ qQ (same generation) (8) 

qQ — > qQ (different generation) (9) 

qQ — > qQ (different generation) (10) 

qq qq (11) 

qq — > QQ (same generation) (12) 

qq — > QQ (different generation) (13) 

qQ — > qQ (same generation) (14) 

qQ — > qQ (different generation), (15) 



with q and Q referring to quarks of different flavours, limited to u-, d-, s-, c- and 6-type (all 
taken as massless) . While the first four processes (with external gluons) were already computed 
in Ref. [3], the eleven four-quark processes are new to this study^. Besides, these four-quark 

processes can be (soft and collinear) infrared divergent, so that gluon bremsstrahlung effects 
ought to be evaluated to obtain a finite cross section at the considered order. In addition, for 
completeness, we have also included the non-divergent 2 — > 3 subprocesses 

qg qqq (16) 

^Notc that in our treatment we identify the jets with the partons from which they originate. 

^Note that gg gg does not appear through 0{a'gaw) nor do qq' — » QQ', qq' — *■ QQ' and qq' — * QQ'. 
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qg qqq (17) 

qg qQQ (same generation) (18) 
qg — >• qQQ (same generation). (19) 

By recalling that at the typical RHIC-Spin energies (e.g., ^/s = 300 and 600 GeV) the 
quark luminosity is much larger than the gluon one, it is clear that are processes with oncom- 
ing quarks that dominate the phenomenology of jet production here. In contrast, at the LHC 
(\/i — 14 TeV), gluon-induced processes are largely dominant, particularly at low Bjorken— x. 
As for what concerns the processes with external gluons, it is worth noticing that no CP- 
violation occurs at tree-level, so that 0{a'ga^J^f) is the first non-trivial order at which parity 
violation is manifest. Regarding four-quark processes, the following should be noted. Parity- 
violating contributions to channels (5)-(15) are induced already at tree-level, through 0(q;|-^). 
Besides, all four-quark channels also exist through the CP- violating 0(q;sQ;ew) [H], although 
subprocesses (9), (10), (13) and (15) only receive Cabibbo-Kobayashi-Maskawa (CKM) sup- 
pressed contributions at this accuracy (i.e., they mainly proceed via CP-conserving 0{ag) 
interactions). Furthermore, notice in general that through (9(agQ;w) there are many more 
diagrams available for channels (1)-(15) than via C(q;s) or indeed 0(q;s«ew) and 0{a'^^^). 
Therefore, in terms of parton luminosity, simple combinatorics and power counting, one should 
expect the impact of 0{aga-w) terms to be large, certainly in parity-violating observables and 
possibly in parity- conserving ones as well. This is what we set out to test in this paper^ for 
the case of RHIC and LHC. (Sec Rcfs. [12, 13] for an account of these effects at Tevatron.) 

Before proceeding further we ought to clarify at this stage that we have only computed 
purely weak effects at one-loop level through 0{ala^^), while in the case of tree-level processes 
via 0{asaEw) and 0{a^y^) also the Electro-Magnetic (EM) contributions are included (and 
so are the interference effects between the two). This is why we are referring in this paper 
to the purely weak terms by adopting the symbol aw, while reserving the notation oew for 
the full EW corrections. Here then, we will have aw = ^em/ sin^ 6'w (with q;em the Electro- 
Magnetic (EM) coupling constant and ^w the weak mixing angle) while a^w will refer to the 
appropriate composition of QED and weak effects as dictated by the SM dynamics. 

We have not computed one-loop EM effects for two reasons. Firstly, their computation 
would be technically very challenging, because the photon in the loop can become infrared (i.e., 
soft and coUinear) divergent, thus requiring also the inclusion of photon bremsstrahlung effects, 
other than of gluon radiation. Secondly, 0{alaEM) terms (in the above spirit, q;em signifies 
here only the contribution of purely EM interactions) would carry no parity-violating effects 
and their contribution to parity-conserving observables would anyway be overwhelmed by the 
well known 0{a^) terms [14] (see also [15, 16]). However, notice that we are not including 
these next-to-leading order (NLO) QCD corrections either, as we are mainly interested in 
parity-violating beam asymmetries. 

Since we are considering weak corrections that may be identified via their induced parity- 
violating effects and since we wish to apply our results to the case of polarised proton beams, it 
is convenient to work in terms of helicity Matrix Elements (MEs). Here, we define the helicity 

•^Subprocesses (16) (19) turn out to be numerically negligible at both machines and whichever the observ- 
able, so that we will not consider them in the remainder. 
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amplitudes by using the formalism discussed in Ref. [17]. At one-loop level such helicity 
amplitudes acquire higher order corrections from: (i) self-energy insertions on the fermions 
and gauge bosons; (ii) vertex corrections and (iii) box diagrams. The expressions for each of 
the corresponding one-loop amplitudes have been calculated using FORM [18] and checked 
by an independent program based on FeynCalc [19]. Internal gauge invariance tests have also 
been performed. The full expressions for the contributions from these graphs are however too 
lengthy to be reproduced here. 

As already mentioned, infrared divergences occur when the virtual or real (bremsstrahlung) 
gluon is either soft or collinear with the emitting parton and these have been dealt with by 
using the formalism of Ref. [20], whereby corresponding dipole terms are subtracted from 
the bremsstrahlung contributions in order to render the phase space integral free of infrared 
divergences. The integrations over the gluon phase space of these dipole terms were performed 
analytically in d-dimensions, yielding pole terms which cancelled explicitly against those of 
the virtual graphs. There remains a divergence from the initial state collinear configuration, 
which is absorbed into the scale dependence of the Parton Distribution Functions (PDFs) and 
must be matched to the scale at which these PDFs are extracted. Recall that the remnant 
initial state collinear divergence at 0{as) is absorbed by the LO dependence of the PDFs. 
Therefore, to 0{alay^), it is sufficient, for the purpose of matching these divergences, to 
consider the LO PDFs. It is also consistent to use the values of the running as obtained form 
the one-loop /9-function. In order to display the corrections due to genuine weak interactions 
the same PDFs and strong coupling are used in the LO and NLO observables. 

The self-energy and vertex correction graphs contain ultraviolet divergences that have 
been subtracted here by using the 'modified' Dimensional Reduction (DR) scheme at the 
scale n ~ Mz- The use of DR, as opposed to the more usual 'modified' Minimal Subtraction 
(MS) scheme, is forced upon us by the fact that the W- and Z-bosons contain axial couplings 
which cannot be consistently treated in ordinary dimensional regularisation. Thus the values 
taken for the running as refer to the DR scheme whereas the EM coupling, aEM, has been 
taken to be 1/128 at the above subtraction point. The one exception to this renormahsation 
scheme has been the case of the self-energy insertions on external fcrmion lines, which have 
been subtracted on mass-shell, so that the external fermion fields create or destroy particle 
states with the correct normalisation 

The top quark entering the loops in reactions with external &'s has been assumed to have 
mass rrit — 175 GeV and width Fj = 1.55 GeV. The Z mass used was Mz — 91.19 GeV and was 
related to the H^-mass, My/, via the SM formula M\y = Mz cos6w, where sm29\Y — 0.232. 
(Corresponding widths were Tz = 2.5 GeV and F^k = 2.08 GeV.) For the strong coupling 
constant, as, we have used the one-loop expression with A^"^"* chosen to match the value 
required by the (LO) PDFs used. The latter were Gehrmann-Stirling set A (GSA) [21] and 
Gliick-Reya-Stratmann- Vogelsang standard set (GRSV-STN) [22]. 

The following beam asymmetries, e.g., can be defined, depending on whether one or both 
beams are polarised: 



All da = 



da. 



++ 




+ 

ALda = 



da 



da^ 



4 



> 



si 



b 
■d 



GSA 

solid/dotted 



jet— production {\ri\ < 1) 
6 = (NLO-LO)/LO 



fi = Et/2 



10^ 



10 
10 

10" 



12 



2.5 
0.0 
-2.5 
-5.0 
-7.5 
-10.0 



-1^ 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 i-i- 


0.5 


-1 1 1 








0.4 




1 1- 






0.3 




1 - 






0.2 




-l 1 1 1 1 1 1 1 


1 1 1 \ 1 1 1 1 1 1 1 1 M 1 1 IM 1 1 


0.1 
0.0 




hJ-1-=FH^'i^ 1 1 fl 1 1 1 1 1 1 1 = 










=11111111 


1111111111111111= 


60 


zl 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 Iz 






40 






^1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 i~l77TT"^rrri^ 


20 



= 1 1 


K^r 1 1 1 1 M>H 1 1 1 1 T 1 1"^ 1 ri--Ki= 




200 250 
T (GeV) 



Figure 1: The dependence of the cross section as well as of the beam asymmetries on the jet 
transverse energy at tree-level (large frames) and the size of the one-loop weak corrections (small 
frames), at the two RHIC-Spin energies ^/s = 300 and 600 GeV. Notice that the pseudo-rapidity 
range of the jets is limited to \r]\ < 1 and the standard jet cone requirement AR > 0.7 is imposed as 
well (although we eventually sum the two- and three-jet contributions). We use GSA as PDFs and 
fi = Et/2 as factorisation/renormalisation scale. Corresponding results for other two energy options 
(mentioned later on) ^/s = 200 and 500 GeV are also given. 



Apyda = da — (i(T_|__|_. (20) 

The first is parity-conserving while the last two are parity- violating^. 

Figs. 1 shows the size of the C(agaw) effects relatively to the well known LO results, 
the latter being defined as the sum of all 0{ag), C(asaEw) and 0(q;|-\y) contributions, for 
the case of RHIC, for two reference energies. Both the differential cross section and the 
above beam asymmetries are plotted, each as a function of the jet transverse energy. The 
(9(agaw) corrections are already very large at cross section level, by reaching — 5(— 9)% at 

= 300(600) GeV, in the vicinity of Et = 120(240) GeV. Effects onto the All asymmetry 
are even larger, with maxima of ~ 25(60)% for Ej- ~ 70(140) GeV, again, in correspondence 
of -s/s = 300(600) GeV. In the case of both Al and Apv, in regions away from the threshold 
at Et ~ Mw/4: (where resonance effects emerge), there is no local maximum for positive or 
negative corrections, as both grow monotonically to the level of +100% (at low Et) and —50 
to —70% (at high Et and with increasing collider energy). All such effects should comfortably 
be observable at RHIC, for the customary values of integrated luminosity, of 200 and 800 pb~^, 
in correspondence of ^/s = 300 and 600 GeV [1]. 

''in the numerical analysis which follows we will assume 100% polarisation of the beams. 
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Figure 2: The dependence of the cross section as well as of the beam asymmetries on the jet 
transverse energy at tree-level (large frames) and the size of the one-loop weak corrections (small 
frames), at the LHC energy ^/s = 14 TeV. Notice that the pseudo-rapidity range of the jets is 
limited to jryj < 2.5 and the standard jet cone requirement AR > 0.7 is imposed as well (although 
we eventually sum the two- and three-jet contributions). We use GSA as PDFs and /j, = Et/2 as 
factorisation / renormalisation scale. 



At the LHC with polarised beams (but standard energy = 14 TeV), the 0(ala^) 
corrections to the total cross section as well as the CP-conserving asymmetry arc reasonably 
under control. In fact, they grow monotonically and reach the ^ —3% and ^ 4% at the 
kinematic limit of the jet transverse energy (as defined by the PDFs), respectively. However, 
it is debatable as to whether these effects can actually be disentangled, as we expect systematic 
experimental uncertainties to be of the same order. Away from the threshold at Ex ~ 
(^(agaw) effects onto the parity-violating asymmetries are instead enormous, as they yield a 
i^'-factor increasing from —2 to —4.5, as Et varies from 80 to 500 GeV. Despite the absolute 
value of the CP-violating asymmetries is rather small in the above interval, the huge LHC 
luminosity (10 fb~^ per year should be feasible for, say, a 70% polarisation per beam [23]) 
would render the above higher order corrections manifest. 

It is intriguing to understand the different behaviours of the 0(aga^j^) effects depending on 
the observable and the collider being considered. To this end, we have presented in Tabs. 1-2 
the contributions to the Et dependent cross section of subprocesses (1)-(15) through 0{alaw) 
separately, at both RHIC-Spin and LHC. The purpose of these tables is to illustrate that 
the leading partonic composition of the 0{agay^) corrections is markedly different at the two 
machines. While at RHIC the key role is played by subprocess (7), at the LHC the conspicuous 
rise of the gluon-luminosity enhances in turn the yield of channel (3) to a level comparable 
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to that of mode (7)^. (The hierarchy among the subprocesses seen in Tabs. 1-2, for fixed jet 
transverse energy, is characteristic across most of the available Ex range at both colliders.) 
The 0{ala^) corrections are particularly large for subprocesses (7)-(8) and (12), mainly in 
virtue of the large combinatorics involved at loop level (as intimated earlier), with respect to 
the LO case. 

Table 1: The contributions of subprocesses (1)-(15) to order agCtw with respect to the fuU LO 
result for the total cross section at RHIC-Spin, at Et = 70 GeV for ^/s = 300 GeV and Et = 140 
GeV for ^/s = 600 GeV. Here, we have paired together the channels with identical Feynman diagram 
topology. We use GSA as PDFs and fi = Et/2 as factorisation/rcnormalisation scale. Column (a) is 
the percentage contribution from the O(agaw) corrections, column (b) is the percentage correction 
to the tree-level partonic subprocess and column (c) is the percentage contribution at the tree-level 
of that partonic subprocess to the differential cross section at the relevant Et- 



Subprocess 


^ = 300 GcV 


^ = 600 GeV 




(a) 


(b) 


(c) 


(a) 


(b) 


(c) 


99 99 






1.35 






1.17 


(1) 


1.73E-05 


0.0267 


0.0650 


-0.000101 


-0.179 


0.0565 


(2) 


4.72E-05 


0.0253 


0.188 


-0.000331 


-0.179 


0.184 


(3)-(4) 


-0.0144 


-0.0589 


24.7 


-0.0608 


-0.264 


23.0 


(5)-(6) 


-0.411 


-0.899 


46.1 


-1.42 


-3.00 


47.1 


(7)-(8) 


-1.57 


-6.64 


23.8 


-3.65 


-14.7 


24.6 


(9)-(10) 


0.000393 


0.0551 


0.718 


-0.00602 


-0.843 


0.711 


(11) 


-0.00236 


-0.252 


0.946 


-0.00664 


-0.714 


0.927 


(12) 


0.00664 


10.90 


0.0613 


0.0146 


25.5 


0.0571 


(13) 


0.00222 


1.23 


0.182 


0.00199 


1.17 


0.169 


(14) 


0.0406 


3.20 


1.28 


0.0282 


2.24 


1.25 


(15) 


0.000350 


0.0491 


0.717 


-0.00572 


-0.803 


0.708 


qq' ^ QQ' or qq' ^ QQ' 






0.00710 






0.0148 


qq' ^ QQ' 






0.00231 






0.00131 


Total 


-1.94 






-5.11 







The different behaviours seen in Figs. 1-2 can easily be interpreted in terms of the LO 
contributions. In this respect, as already mentioned. Tabs. 1-2 clearly make the point that the 
jet phenomenology at RHIC-Spin is dominated by subprocesses initiated by quarks only while 
at the LHC gluons- induced channels arc generally predominant. At RHIC energies, LO pro- 
duction through order is dominated by channel (5) whereas at the LHC the overwhelmingly 
dominant ag channels are gg — > gg (which is not subject to C(agQ;w) corrections, as already 
mentioned) and subprocess (3). Besides, channel (5) through ^(ag) is mainly concentrated 

^The relevance of the latter throughout originates from the combination of a always sizable valence quark 
luminosity and a large Feynman diagram combinatorics, as opposed to, e.g., a gluon luminosity steeply in- 
creasing with the collider energy but combined with a small numbers of graphs [3]. 
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Table 2: The contributions of subprocesses (1)-(15) to order agaw with respect to the full LO 
result for the total cross section at LHC, at Et = 300 GeV for ^/s = 14 TeV. Here, we have 
paired together the channels with identical Feynman diagram topology. We use GSA as PDFs and 
fi = Et/2 as factorisation/renormalisation scale. Column (a) is the percentage contribution from the 
0{aga\y) corrections, column (b) is the percentage correction to the trcc-lcvel partonic subprocess 
and and column (c) is the percentage contribution at the tree-level of that partonic subprocess to 
the differential cross section at the relevant Et- 



Subprocess 


v/i = 14 TeV 




(a) 


(b) 


(c) 








41.9 


(1) 


-0.0315 


-1.25 


1.89 


(2) 


-0.00386 


-1.27 


0.228 


(3)-(4) 


-0.455 


-0.711 


47.8 


(5)-(6) 


-0.112 


-4.92 


1.70 


(7)-(8) 


-0.431 


-17.3 


1.87 


(9)-(10) 


-0.0330 


-2.67 


0.926 


(11) 


-0.0328 


-1.85 


1.328 


(12) 


0.0466 


64.8 


0.0540 


(13) 


0.0031(3 


1.50 


0.158 


(14) 


0.0131 


0.821 


1.196 


(15) 


-0.0325 


-2.64 


0.924 


qq' ^ QQ' or qq' ^ QQ' 






0.0221 


qq' ^ QQ' 






0.000979 


Total 


-1.075 







at low Et while with growing Et the C(q;sq;ew) and - particularly - 0{a^^^) terms gain 
in relative importance. Furthermore, 0(«s) terms entering channel (5) do not contribute, 
obviously, to the parity-violating asymmetries. Therefore, it should not be surprising to see 
at RHIC-Spin that our corrections are very large in the case of the latter, where the LO term 
is 0{asOi-Ew), respect to which the corrections computed here are suppressed only by one 
power of as. We attribute instead the size of the 0{ala^^) effects on the cross section and the 
parity-conserving asymmetry again to the fact that through (9(aga-vv) there are many more 
diagrams available for such channels than via ^(0:3) or indeed 0{asC(Ew) and (9(q;|-^). As 
for the LHC, the fact that gg — > gg and subprocess (3) vastly dominates through 0{al) the 
da/dET distribution explains why 0{agay/) corrections are limited to the percent level. In 
All, which has no gg — > gg component, 0{aga-w) effects become somewhat more visible in 
comparison. 

Furthermore, in the case of the LHC, one should note the monotonic rise of the corrections 
with increasing jet transverse energy, for all observables studied, which can be attributed to 
the so-called Sudakov (leading) logarithms [24, 25] of the form ctw log^ (ii^|./M^), which appear 
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in the presence of higher order weak corrections. These 'double logarithms' are due to a lack 
of cancellation of infrared (both soft and coUinear) virtual and real emission in higher order 
contributions due to 1^-exchange, arising from a violation of the Bloch-Nordsieck theorem 
occurring in non-Abelian theories. (In fact, if events with real Z radiation are vetoed in 
the jet sample, aw log^(i?T/^I) terms would also affect the corrections [12].) Clearly, at 
LHC energies, Et can be very large, thus probing the kinematic regime of these logarithmic 
effects, which instead affected RHIC only very mildly. Combine then the effects of such large 
logarithms with the fact that Ai^ and Apy receive no pure QCD contributions, and one can 
explain the enormous (and increasing with Et) Oia^ayi) corrections to these two obscrvables. 
In fact, recall that another way of viewing the (9(agQ;w) terms computed here is as first order 
QCD corrections to the ^(asttw) terms, which are the leading order contributions to Ai^ and 
Apy. Prom this perspective then, the large results reported here can be understood as large 
0(q;s) corrections. Furthermore, also recall here the following two aspects, already mentioned. 
Firstly, there are several partonic processes which are CKM suppressed at C(q;sq;ew) but 
which occur without CKM suppression at (9(aga;w)- Secondly, processes involving external 
gluons and weak interactions occur for the first time at ^(agaw)- 

As one of the purposes of polarised colliders is to measure polarised structure functions, 
in the ultimate attempt to reconstruct the proton spin, it is of some relevance to see how 
the Oipcyxyq) results obtained so far for GSA compare against GRSV-STN. This is done 
in Figs. 3-4, where we have also adopted the different choice ^ = EauiEr) as factorisa- 
tion/renormalisation scale, i.e., the centre-of-mass energy at parton level v^(jet transverse 
energy), for the GSA(GRSV-STN) set. A comparison between the GSA curves in Figs. 3-4 
and those in Figs. 1-2 reveals that the scale dependence of our corrections is not very substan- 
tial for a given PDF set (the same is true for the case of GRSV-STN). In contrast, depending 
on the choice of PDFs, corrections through Oia^a-^) can be very different for each obscrvables 
studied at both RHIC-Spin and LHC, with the exception of the cross section in cither case. 

Altogether, the results presented here point to the extreme relevance of one-loop Oia^^a^) 
weak contributions for precision analyses of jet data produced in polarised proton-proton 
scattering at RHIC. We have confirmed that this would be the case also at a polarised LHC, 
which has been discussed as one of the possible upgrades of the CERN collider. The size of 
the afore-mentioned corrections, relative to the lowest order results, is rather insensitive to 
the choice of factorisation/renormahsation scale, yet it shows some sizable dependence on the 
polarised PDFs used. EM effects were neglected here because they are not subject to parity- 
violating effects. However, their computation is currently in progress. The inclusion of NLO 
terms from pure QCD, through (^(ag), is also in order, as they can produce effects of order 
100%, even to the parity-violating asymmetries, though in this case they will only amount to 
a rescaling (within a factor of 2 at the most) of the normalisation, not to a change in shape. 
We are now working towards the full 0{p^ results including beam polarisation effects [26]. 

Finally, extrapolation of our results to other collider energies, chiefly for RHIC, for oper- 
ation at = 200 and 500 GeV, is straightforward. As expected, these results do not differ 
dramatically from those obtained for 300 and 600 GeV, respectively. Fig. 1 also shows the size 
of the corrections at these two additional energies, for our usual obscrvables and default choice 
of PDFs and factorisation/renormahsation scale. For the purpose of emulating the effects of 
0(q;sq;^) terms at whichever collider and energy, we make available our code on request. 
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Figure 3: The dependence of the corrections to the cross section as weU as the beam asymmetries on 
the jet transverse energy for two sets of PDFs, GSA and GRSV-STN, at the two RHIC-Spin energies 
^/s = 300 (curves extending to 150 GeV) and 600 GeV (curves extending to 300 GeV). Notice that 
the pseudo-rapidity range of the jets is limited to |??| < 1 and the standard jet cone requirement 
Ai? > 0.7 is imposed as weh (although we eventually sum the two- and three-jet contributions). We 
use /X = EcmiEr) as factorisation/renormalisation scale in the case of GSA(GRSV-STN). 
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